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Abstract 

In general, the behavior of large and complex aggregates of elementary 
components can not be understood nor extrapolated from the properties of 
a few components. The brain is a good example of this type of networked 
systems where some patterns of behavior are observed independently of 
the topology and of the number of coupled units. Following this insight, we 
have studied the dynamics of different aggregates of logistic maps accord- 
ing to a particular symbiotic coupling scheme that imitates the neuronal 
excitation coupling. All these aggregates show some common dynami- 
cal properties, concretely a bistable behavior that is reported here with 
a certain detail. Thus, the qualitative relationship with neural systems 
is suggested through a naive model of many of such networked logistic 
maps whose behavior mimics the waking-sleeping bistability displayed by 
brain systems. Due to its relevance, some regions of multistability are 
determined and sketched for all these logistic models. 
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1 INTRODUCTION 



The brain is a natural networked system [1] [2] . The understanding of this com- 
plex system is one of the most fascinating scientific tasks today, concretely how 
this set of millions of neurons can symbiotically interact among them to give 
rise to the collective phenomenon of human thinking [3], or, in a simpler and 
more realistic approach, what neural features can make possible, for example, 
the birdsongs [I]. Different aspects of neurocomputation take contact on this 
problem: how brain stores information and how brain processes it to take de- 
cisions or to create new information. Other universal properties of this system 
are more evident. One of them is the existence of a regular daily behavior: the 
sleep- wake cycle [5j |6] . 

NIGHT 
1 

SLEEP 
state 

Figure 1: Brain bistability provoked by the solar light cycle. 

The internal circadian rhythm is closely synchronized with the cycle of sun 
light (Fig. [T]). Roughly speaking and depending on the particular species, the 
brain is awake during the day and it is slept during the night, or vice versa. All 
mammals and birds sleep. There is not a well established law relating the size of 
the animal with the daily time it spends sleeping, but, in general, large animals 
tend to sleep less than small animals (Fig. 2). Hence, at first sight, the emergent 
bistable sleep- wake behavior seems not depend on the precise architecture of the 
brain nor on its size. This structural property would mean that, if we represent 
the brain as a complex assembly of units HI [9] , the possible bistability, where 
large groups of neurons can show some kind of synchronization, should not 
depend on the topology (structure) nor on the number of nodes (size) of the 
network (Fig. |3j). 



ANIMAL 
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WEIGHT 

(gr.) 


SLEEP 
TIME 

(hours) 


Cat 


30 


13 


Dog 


70 


11 


Human (-) 


350 


16 


Human (+) 


1300 


8 


Horse 


530 


2 


Elephant 


4500 


4 



Figure 2: Brain size and mean sleeping time for different animals. (Humans(-) 
represent new born humans and Humans (+) represent middle age humans). 

Then, it is essential the type of local dynamics and the excitation/inhibition 
coupling among the nodes that must be implemented in order to get bistability 
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The waking and sleeping brain 
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Figure 3: The waking-sleeping brain. 



as a possible dynamical state in a complex network. So, on one side, it has been 
recently argued [10] that the distribution of functional connections p(k) in the 
human brain, where p(k) represents the probability of finding an element with k 
connections to other elements of the network, follows the same distribution of a 
scale- free network [llj . that is a power law behavior, p(k) ~ k^ 1 , with 7 around 
2. This finding means that there are regions in the brain that participate in 
a large number of tasks while most of the regions are only involved in a tiny 
fraction of the brain's activities. On the other side, it has been shown by 
Kuhn et al. [12] the nonlinear processing of synaptic inputs in cortical neurons. 
They studied the response of a model neuron with a simultaneous increase of 
excitation and inhibition. They found that the firing rate of the model neuron 
first increases, reaches a maximum, and then decreases at higher input rates. 
Functionally, this means that the firing rate, commonly assumed to be the 
carrier of information in the brain, is a non-monotonic function of balanced 
input. These findings do not depend on details of the model and, hence, they 
are relevant to cells of other cortical areas as well. 



COMPLEX 
NETWORK 



NONLINEAR 
OSCILLATORS 



BISTABILITY 



YES 



Figure 4: QUESTION: Is it possible to implement some kind of coupling and 
nonlinear dynamics in each node of a complex network in order to get bistabil- 
ity?. ANSWER: Yes. 



Putting together all these facts, we arrive to the central question that we 
want to bring to the reader: Is it possible to reproduce the bistability in a 
complex network independently of the topology and of the number of nodes?. 
The answer is 'Yes' (Fig. rjj. What kind of local dynamics and coupling among 
nodes must be implemented in order to get this behavior?. In Section [21 we give 
different possible strategies for the coupling and the local dynamics which should 
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be implemented in a few coupled functional units in order to get a bistable 
behavior. Then, in Section [3j the same model is exported on a many units 
network where the desired bistability is also retained. In view of the results, the 
possibility of constructivism in the world of complex systems in general, and in 
the neural networks in particular, is suggested in the Conclusion. 

2 MODELS OF A FEW COUPLED 
FUNCTIONAL UNITS 

2.1 General model 

Our approach considers the so called functional unit, i.e. a neuron or group 
of neurons (voxels), as a discrete nonlinear oscillator with two possible states: 
active (meaning one type of activity) or not (meaning other type of activity). 
Hence, in this naive vision of the brain as a networked system, if x l n , with 
< x l n < 1, represents a measurement of the ith functional unit activity at 
time n, it can be reasonable to take the most elemental local nonlinearity, for 
instance, a logistic evolution |13) . which presents a quadratic term, as a first 
toy-model for the local neuronal activity: 

x n+l = Pi x% n\\ ~ x n)- W 



Figure 5: Discrete nonlinear model for the local evolution of a functional unit. 

It presents only one stable state for each pi . Then, there is no bistability in 
the basic component of our models. For < pi < 1, the dynamics dissipates 
to zero, x l n = 0, then it can represent the functional unit with no activity. For 
1 < pi < 4, the dynamics is non null and it would represent an active functional 
unit. 

We can suppose that this local parameter pi is controlled by the signals 
of neighbor units, simulating in some way the effect of the synapses among 
neurons. Excitatory and inhibitory synaptic couplings have been shown to be 
determinant on the synchronization of neuronal firing. For instance, facilitatory 
connections are important to explain the neural mechanisms that make possible 
the object representation by synchronization in the visual cortex [14] . While 
it seems clear that excitatory (symbiotic) coupling can lead to synchronization, 
frequently inhibition rather than excitation synchronizes firing [TS] . The impor- 
tance of these two kinds of coupling mechanisms has also been studied for other 
types of neurons, v.g., motor neurons [16] . 

If a neuron unit simultaneously processes a plurality of binary input signals, 
we can think that this local information processing is reflected by the parameter 
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pi. The functional dependence of this local coupling on the neighbor states is 
essential in order to get a good brain-like behavior (i.e., as far as the bistability 
of the sleep-wake cycle is concerned) of the network. As a first approach, we 
can take pi as a linear function depending on the actual mean value, X. % n , of the 
neighboring signal activity and expanding the interval (0, 4) in the form: 

Pi — Pi (3^Q + 1), (excitation coupling) (2) 
or 

pi — pi (—3X^ + 4), (inhibition coupling) (3) 

with 

Ni 

K = ^E<- ( 4 ) 

Ni is the number of neighbors of the ith functional unit, and pi, which gives us 
an idea of the interaction of the functional unit with its first-neighbor functional 
units, is the control parameter. This parameter runs in the range < pi < p m ax, 
where p ma x > 1- When pi = p for all i, the dynamical behavior of these 
networks with the excitation type coupling [3 [8] presents an attractive global 
null configuration that has been identified as the turned off state of the network. 
Also they show a completely synchronized non-null stable configuration that 
represents the turned on state of the network. Moreover, a robust bistability 
between these two perfect synchronized states is found in that particular model. 
For different models with a few coupled functional units we sketch in the next 
subsections the regions where they present a bistable behavior. The details of 
the complete unfolding 17 of these dynamical systems can be found in the 
references [HHni [20] ■ 

2.2 Models of two functional units 

Let us start with the simplest case of two interconnected (x n ,y n ) functional 
units. Three different combinations of couplings are possible: (excitation, excitation) , 
(excitation, inhibition) and (inhibition, inhibition). 

I 1 excitation i 1 

1 1 inhibition I 

Figure 6: Two functional coupled units. 

As we are concerned with interactions with a certain degree of symbiosis 
between the coupled units, the regions of parameter space where bistability is 
found for the first two cases are presented here. Experimental systems where 
similar couplings have been found or implemented for two cells systems have 
been reported in the literature [2TJ [22| ■ 
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2.2.1 Model with mutual excitation 

The dynamics of the (excitation, excitation) case |18l is given by the coupled 
equations: 

x„+i = p (3y„ + l)x n (l - x n ), (5) 
Vn+i = P (3a;„ + l)y«(l — y n )- (6) 
The regions of the parameter space (Fig. [7]) where we can find bistability are: 

• For 0.75 < p < 0.86, the synchronized state, x+ = (x,x) = Pi, with 
a: = i{l + (4— 1)^}, which arises from a saddle- node bifurcation for the 
critical value p = 0.75, is a stable turned on state. This state coexists 
with the turned off state xg = 0. The system presents now bistability 
and depending on the initial conditions, the final state can be xg or x+. 
Switching on the system from xg requires a level of noise in both functional 
units sufficient to render the activity on the basin of attraction of 2+ . On 
the contrary, switching off the two functional units network can be done, 
for instance, by making zero the activity of one functional unit, or by 
doing the coupling p lower than 0.75. 

• For 0.86 < p < 0.95, the active state of the network is now a period-2 oscil- 
lation, namely the period-2 cycle (Ps,Pq) in Fig. [7l This new dynamical 
state bifurcates from x+ for p — 0.86. A smaller noise is necessary to 
activate the system from xg. Making zero the activity of one functional 
unit continues to be a good strategy to turn off the network. 

• For 0.95 < p < 1, the active state acquires a new frequency and presents 
quasiperiodicity (the invariant closed curves of Fig. [7])- It is still possible 
to switch off the network by putting to zero one of the functional units. 



1 = fl.a A= 0.9 3. = 9.9S 




Figure 7: Bistability in 2 functional units with excitation type coupling (p = A) . 
2.2.2 Model with excitation + inhibition 

The dynamics of the (excitation, inhibition) case |19) is given by the coupled 
equations: 

x n +i = p (3y„ + l)x„(l - x n ), (7) 
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y n+ i = p (-3x„ + 4)j/ n (l-t/„). 



(8) 



The regions of the parameter space (Fig. [5]) where we can find bistability are: 

• For f.051 < p < 1.085f, a stable period three cycle (Qi,Q2,Qs) appears 
in the system. It coexists with the fixed point P4. When p is increased, 
a period-doubling cascade takes place and generates successive cycles of 
higher periods 3-2". The system presents bistability. Depending on the 
initial conditions, both populations (x n ,y n ) oscillate in a periodic orbit 
or, alternatively, settle down in the fixed point. The borders between the 
two basins are complex. 

• For 1. 0851 < p < 1.0997, an aperiodic dynamics is possible. The period- 
doubling cascade has finally given birth to an order three cyclic chaotic 
band(s) (^-31,^32,^33). The system can now present an irregular oscilla- 
tion besides the stable equilibrium with final fixed populations. The two 
basins are now fractal. 




Figure 8: Bistability in 2 functional units with excitation+inhibition type cou- 
pling (p = A). 



2.3 Models of three functional units 

Following the strategy given by relation ([2|[3]) several models with three func- 
tional units can be established. We have studied in some detail three of them 
[20] and their bistable behavior is reported here. 

2.3.1 Model with local mutual excitation 

Let us start with the case of three alternatively interconnected (x n , y n , z n ) func- 
tional units under a mutual excitation scheme. 

Then the dynamics of the system is given by the coupled equations: 

x n+ i = p (3y„ + l)x n (l - x n ), (9) 
y n+1 = p(3z n + l)y n (l-y n ), (10) 
Zn+i = P (3a;„ + l)z n (l - z„). (II) 
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Figure 9: Three alternatively coupled functional units under the excitation 
scheme. 





Figure 10: Bistability in 3 functional units with local excitation type coupling 
(p = A). 

The regions of the parameter space where we have found bistability are: 

• For 0.93310 < p < 0.95334, a big invariant closed curve (ICC) CI coex- 
ists with a period-3 orbit that bifurcates, first to an order 3-cyclic ICC 
(Fig. [TO)) , and finally to an order-3 weakly chaotic ring (WCR) before 
disappearing. 

• For 0.98418 < p < 0.98763, the ICC CI coexists with another ICC C2 (see 
Ref. 20 ) that becomes chaotic, by following a period doubling cascade 
of tori, before disappearing. 

• For 1.00360 < p < 1.00402, the ICC CI coexists with a high period 
orbit that gives rise to an ICC C3. This ICC also becomes a chaotic 
band (see Ref. [20]) by following a period doubling cascade of tori before 
disappearing. 

2.3.2 Model with global mutual excitation 

We expose now the case of three globally interconnected (x n , y n , z n ) functional 
units under a mutual excitation scheme. 

Then the dynamics of the system is given by the coupled equations: 

X n +1 = p (x n + y n + Z n + l)x n (l - X n ), (12) 
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Vn+1 = P {X n + Vn + Z n + l)j/„(l - y n ), (13) 
Zn+1 = P (x n + Vn + Z n + l)z„(l - Z n ) . (14) 



| | global mutual 

^ ^ ^ 

r—K A ^-X excitation 



Figure 11: Three globally coupled functional units under the excitation scheme. 




Figure 12: Bistability in 3 functional units with global excitation type coupling 
(p = A). 

For the whole range of the parameter, < p < 1.17, bistability is present in 
this system: 

• Firstly, two order-2 cyclic ICC coexist before becoming two order-2 cyclic 
chaotic attractors (Fig. [T2|) by contact bifurcations of heteroclinic type. 
Finally the two chaotic attractors become a single one before disappearing. 

2.3.3 Model with partial mutual excitation 

The new case [5] of three partially interconnected (x n , y n , z n ) functional units 
under a mutual excitation scheme is represented in Fig. 1131 

The dynamics of the system is given by the coupled equations: 

x„+i = p (3(y» + z n )/2 + l)x„(l - x n ), (15) 
Vn+i = p (3(xn + z n )/2 + l)j/„(l - y n ), (16) 

Zn+l = P (3(Xn + y n )/2 + l)z n (l - Z n ). (17) 
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EMJial mutual 
excitation 



Figure 13: Three partially coupled functional units under the excitation scheme. 



a=1 .05 




Figure 14: Multistability in 3 functional units with partial excitation type cou- 
pling. The system presents three order-2 ICC for p = a = 1.05. 



The rough inspection of this system puts in evidence the existence of different 
regions of multistability in the parameter space. These are: 

• For 0.93 < p < 1.04, there is coexistence among three cycles of period-2. 

• For 1.04 < p < 1.06, the three cycles bifurcate giving rise to three order-2 
ICC (Fig. m. 

• For 1.06 < p < 1.08, the system can present three mode-locked periodic 
orbits, each one with period multiple of 6, displaying period doubling 
cascades giving rise to three chaotic cyclic attractors. Three chaotic cyclic 
attractors of order 2 are also possible in this region (Fig. fTS"]) . 

• For p > 1.08, the chaotic cyclic attractors collapse in an unique chaotic 
attractor (Fig. IT6]) . 

Also, other multistable situations can be found for some particular values of the 
parameter p in the former intervals, such as it can be seen in Fig. [T71 where the 
^-projection of a generic bifurcation diagram is plotted. Similar diagrams are 
found for the y- and ^-projections. 
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1=1.078 




Figure 15: Multistability of three chaotic cyclic attractors of order 2 for p — 
a = 1.078. 

3 MODEL OF MANY COUPLED 
FUNCTIONAL UNITS 

3.1 The model 

The complete synchronization [531 123] of the network means that x l n — x n for 
all i , with i = 1,2, ... ,N and N J> 1 25, 26 . In this regime, we also have 
X n — x n . The time evolution of the network [8] on the synchronization manifold 
is then given by the cubic mapping: 

x n+ i = p (3x n + 1) x„(l - Xn). (18) 

The fixed points of this system are found by solving x„+i = x n . The solutions 
are xg = and x± — ^{1 ± (4— The first state xg is stable for < p < 1 

and x± take birth after a saddle-node bifurcation for p = po = 0.75. The node 
x+ is stable for 0.75 < p < 1.157 and the saddle x_ is unstable. Therefore 
bistability between the states 

x n =x e , Vi — > TURNED OFF STATE, (19) 
x i n =x+, Vi — > TURNED ON STATE, (20) 

seems to be also possible for p > po — 0.75 in the case of many interacting units. 
But stability on the synchronization manifold does not imply the global stability 
of it. Small transverse perturbations to this manifold can make unstable the 
synchronized states. Let us suppose then a general local perturbation Sx l n of 
the element activity, 

x n = x* + 8x n , (21) 
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Figure 16: An unique chaotic attractor for p — a = 1.09. 



with x* representing a synchronized state, xg or x+. We define the perturbation 
of the local mean-field as 

SXi = w ^2Sxl (22) 

1 3 = 1 

If these expressions are introduced into equation ([1]) , the time evolution of the 
local perturbations are found: 

Sx l n+1 = p(3x» + 1)(1 - 2x*)6x i n + px*(l - x*)6Xi. (23) 

The dynamics for the local mean-field perturbation is derived by substituting 
this last expression in relation (|22"j) . We obtain: 

1 N * 

5X* n+1 = p (3x* + 1)(1 - 2x*)5Xl l + 3px*(l - x*) — £ 5X1 (24) 

We express now the local mean- field perturbations of the first- neighbors as func- 
tion of the local mean-field perturbation SX^ by defining the local operational 
quantity cr™, 

N t 

— ^6X1 = ^6X1 (25) 

' 3=1 

which is determined by the dynamics itself. If we put together the equations 
l]23H24p . the linear stability of the synchronized states holds as follows: 

/ 5x l n+1 \ _f p(3x* + 1)(1 - 2x») px*(l -x*) \ / Sx l n 

\ 5X l n+1 ) ~ \ p{3x, + 1)(1 - 2a;,) + 3p<T* n x m (l -sc.) J \ SX % n 

(26) ' 

Let us observe that the only dependency on the network topology is included in 
the quantity u l n . The rest of the stability matrix is the same for all the nodes 
and therefore it is independent of the local and global network organization. 
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Model with partial mutual excitation : x versus p 
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Figure 17: Bifurcation diagram projected in the x coordinate for the initial 
conditions (xq — 0.37, yo — 0.36, zq = 0.33) as a function of the parameter p. 



The turned off state is x* = xg = 0. The eigenvalues of the stability matrix 
are in this case Ai = A2 = p. Then, this state is an attractive state in the 
interval < p < 1. It loses stability for p = 1, then the highest value p/ of the 
parameter p where bistability is still possible satisfies pf < 1. 

The turned on state x + verifies a;* = x+ = |{1 + (4 — |) 3}. If we suppose 

a l n = a, the eigenvalues of the stability matrix are Ai = 2 — 2p — p(4 — |)i 

and A2 = Ai + ^(3 — 2p + p(4 — |)')- Let us observe that Ai = — 1 for p = 1. 
This implies that the parameter p c for which the synchronized state x+ looses 
stability verifies p c < 1. Depending on the sign of tr, we can distinguish two 
cases in the behavior of p c : 

• If < a < 1, we find that | A2 |< 1. Then x + bifurcates through a 
global flip bifurcation for p = p c = 1. In this case, the bifurcation of the 
synchronized state x+ for p c = 1 coincides with the loss of the network 
bistability for pf = 1. Hence p c = pf = 1 for this kind of networks, 
and the bistability holds between xg and x+ in the parameter interval 
Po = 0.75 < p < p c = pf = 1. As an example, an all-to-all network shows 
this behavior because a = 1 . This is represented in the inset of Fig. [151 

• If— 1 < a < 0, then A2 = — 1 is obtained for a p — p c smaller than 1. 
Therefore it is now possible to obtain an active state different from x+ in 
the interval p c < p < pf . For instance, simulations show that the global 
flip bifurcation of the synchronized state for a scale free network occurs for 
Pc = 0.87±0.01. A value of p = 0.866 is obtained from the stability matrix 
by taking a = — 1. For this particular network it is also found that pt = 1. 
Then, bistability is possible in the range po = 0.75 < p < pf — 1 for this 
kind of configuration. But now an active state with different dynamical 
regimes is observed in the interval p c = 0.87 < p < p/ = 1. If we identify 
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Figure 18: Stable states (xg,x+) of the network for < p < 1. Let us observe 
the two zones of bistability: Po < p < p c and p c < p < Pf. The main figure 
corresponds to a scale free network made up of TV = 10 elements: po = 0.75, 
Pc = 0.87±0.01 and pf — 1. The inset shows the same graph but in an all-to-all 
network of the same size: po = 0.75, p c — pf — 1. Initial conditions for the x^s 
were drawn from a uniform probability distribution in the interval (0, 1). 



the capacity of information storing with the possibility of the system to 
access to complex dynamical states, then, we could assert, in this sense, 
that a scale free network has the possibility of storing more elaborated 
information in the bistable region that an all-to-all network. 

Let us note that a also indicates a different behavior of local dissipation, as 
expression (|25l) suggests. A positive a means a local in-phase oscillation of the 
node signal and mean-field perturbations. A negative tr is meaning a local out of 
phase oscillation between those signal perturbations. Hence, a also brings some 
kind of structural network information. In all the cases the stability loss of the 
completely synchronized state is mediated by a global flip bifurcation. The new 
dynamical state arising from that active state for p = p c is a periodic pattern 
with a local period-2 oscillation. The increasing of the coupling parameter 
monitors other global bifurcations that can lead the system towards a pattern 
of local chaotic oscillations. 

3.2 Transition between On-Off States 

We proceed now to show the different strategies for switching on and off a 
random scale free network. The choice of this network is suggested by the 
recent work |10j on the connections distribution among functional units in the 
brain. They find it to be a power- law distribution. Following this insight [5], a 
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scale-free network following the Barabasi-Albert (BA) recipe [TT] is generated. 
In this model, starting from a set of mo nodes, one preferentially attaches each 
time step a newly introduced node to m older nodes. The procedure is repeated 
N — tuq times and a network of size N with a power law degree distribution 
P(k) ~ fc -7 with 7 = 3 and average connectivity (A;) = 2m builds up. This 
network is a clear example of a highly heterogenous network, in that the degree 
distribution has unbounded fluctuations when N — > oo. The exponent reported 
for the brain functional network has 7 < 3. However, studies of percolation and 
epidemic spreading [57J [25] on top of scale-free networks have shown that the 
results obtained for 7 = 3 are consistent with those corresponding to lower values 
of 7, with 7 > 2. As explained before, network bistability between the active 
and non active states is here possible in the interval po = 0.75 < p < Pf = 1 

(Fig. ED. 

Switching off the network 

Two different strategies [5J can be followed to carry the network from the active 
state to that with no activity (Fig. [T5)) . 

• Route I: By doing the coupling p lower than p Q . This is the easiest and 
more natural way of performing such an operation. In our naive picture of 
a brain-like system, it could represent the decrease (or increase, it depends 
on the specific function) of the synaptic substances that provokes the 
transition from the awake to the sleep state. The flux of these chemical 
activators is controlled by the internal circadian clock, which is present in 
all animals, and which seems to be the result of living during millions of 
years under the day/night cycle. 

• Route II: By switching off a critical fraction of functional units for a fixed 
p. Evidently, at first sight, this strategy has no relation with the behavior 
of a real brain-like system. Thus, this is done by looking over all the 
elements of the network, and considering that the element activity is set 
to zero with probability A (which implies that on average XN elements are 
reset to zero). The result of this operation is shown in Fig. [TSJ Here, it 
is plotted for different p's the relative size of the biggest (giant) cluster of 
connected active nodes in the network versus A. Note that this procedure 
does not take into account the existence of connectivity classes, but all 
nodes are equally treated. The procedure is thus equivalent to simulations 
of random failure in percolation studies |27) . The strategy in which highly 
connected functional units are first put to zero is more aggressive and 
leads to quite different results. 

Each curve presents three different zones depending on A: 



- the robust phase: For small A, the network is stable and only those 
states put to zero have no activity. There is a linear dependence on 
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Figure 19: Turning off a scale free network. Three different phases in the 
behavior of the giant cluster size versus A (fraction of switched off nodes) are 
observed. These three phases are illustrated for p = 0.96: the robust phase, the 
weak phase and the catastrophic phase (see the text). Other network parameters 
are as those of Fig. [18] 

the giant cluster size with A. In this stage, the switched off nodes 
do not have the capacity to transmit its actual state to its active 
neighbors. 

- the weak phase: For an intermediate A, the nodes with null activity 
can influence its neighborhood and switch off some of them. The 
linearity between the size of the giant cluster and A shows a higher 
absolute value of the slope than in the robust zone. 

- the catastrophic phase: When a critical A c is reached, the system 
undergoes a crisis. The sudden drop in this zone means that a small 
increase of the non active nodes leads the system to a catastrophe; 
that is, the null activity is propagated through all the network and 
it becomes completely down. 

It is worth noticing that when the system is outside the bistability region 
for p > 1, the catastrophic phase does not take place. Instead, the turned off 
nodes do not spread its dynamical state and the neighboring nodes do not die 
out. This is because the dynamics of an isolated node is self-sustained when 
p > 1. Consequently, it is observed that the network breaks down in many 
small clusters and the transition resembles that of percolation in scale free nets 

nana]. 



16 



Switching on the network 

Two equivalent strategies [5] can be followed for the case of turning on the 
network (Fig. 120]) : 

• (I) For a fixed p, we can increase the maximum value e of the noisy signal, 
which is randomly distributed in the interval (0, e) over the whole system. 
When e attains a critical value e c , the noisy configuration can leave the 
basin of attraction of x$, whose boundary seems to have the form in phase 
space of a "hollow cane" around it, and then the network rapidly evolves 
toward the turned on state; 

• (II) If this operation is executed by letting e to be fixed and by increasing 
the coupling parameter p, the final result of switching on the network 
is reached when p takes the value for which e = e c . The final result is 
identical in both cases. 

Let us remark that the strategy equivalent to the former Route II, that is, 
the switching on of a critical fraction of functional units, is not possible in this 
case. It is a consequence of the fact that a switched off functional unit can not 
be excited by its neighbors and it will maintain indefinitely the same dynamical 
state (xi — 0). 

Finally, let us note that, from Fig. 1201 a bigger p requires a smaller e c 
to switch on the net. Observe that this behavior could be interpreted in our 
approach as the smaller level of noise that is needed for awaking a brain-like 
system that it is departing from the sleeping state. 

4 CONCLUSION 

One of the more challenging problems in nonlinear science is the goal of under- 
standing the properties of neuronal circuits |30j . Synchrony and multistability 
are two important dynamical behaviors found in those circuits [31( I32( I33j . In 
this work, different coupling schemes for networks with local logistic dynamics 
are proposed [7]. It is observed that these types of couplings generate a global 
bistability between two different dynamical states [9] . This property seems 
to be topology and size independent. This is a direct consequence of the local 
mean-field multiplicative coupling among the first-neighbors. If a formal and 
naive relationship is established between these two states and the sleep-wake 
states of a brain, respectively, one would be tempted to assert that these types 
of couplings in a network, regardless of its simplicity, give us a good qualitative 
model for explaining that specific bistability. Following this insight, several low- 
dimensional systems with logistic components coupled under these schemes have 
been presented and the regions where the dynamics shows bistability have been 
identified. The extension of this type of coupling to a general network with local 
logistic dynamics has also been achieved. This system presents global bistability 
between an active synchronized state and another synchronized state with no 
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Figure 20: Turning on a scale free network. For a fixed p, a noisy signal randomly 
distributed in the interval (0, e) is assigned to every node. When e reaches the 
critical level e c the network becomes switched on. Other network parameters 
are as those of Fig. [18] 

activity. This property is topology and size independent. This is a direct conse- 
quence of the local mean-field multiplicative coupling among the first-neighbors. 
If a formal and naive relationship is established between the switched off and 
switched on states of that network, and the sleep-wake states of a brain, re- 
spectively, one would be tempted to assert that this model, regardless of its 
simplicity, is a good qualitative representation for explaining that specific bista- 
bility. Furthermore, on more theoretical grounds, the results obtained here 
point out the entangled interplay between topology and function in networked 
systems [34] where complex structures coexist with nonlinear dynamics. 
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